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Hybrid  Electric  Vehicles  (HEVs)  are  becoming  more  popular  than  pure  electric  ones,  nowadays.  This  is 
because  of  their  better  performance,  economic  advantages  and  higher  operating  range.  However,  their 
potential  advantages  extremely  depend  on  their  system  design,  most  importantly  their  battery  system 
design.  Batteries’  life  requirements  as  well  as  the  cost  of  replacing  them  at  the  end  of  their  life  period, 
currently  limit  manufacturers  to  bring  HEVs  into  play,  even  though  their  fuel  economy  reduces  their 
everyday  cost  considerably.  Generally,  inappropriate  discharge/charge  patterns  would  result  in  loss  in 
batteries’  life.  In  the  present  study,  an  optimization  based  control  strategy  has  been  proposed  for  the 
series  HEVs  in  order  to  maximize  the  efficiency  of  the  power-train  while  minimizing  the  loss.  A  genetic 
algorithm  is  implemented  to  optimally  evaluate  the  control  algorithm’s  parameters.  The  approach  is 
then  compared  to  two  main  control  strategies,  namely  thermostatic  control  strategy  and  power  follower 
control  strategy.  The  computational  procedure  of  the  genetic  algorithm  is  discussed,  and  a  simulation 
study  based  on  a  model  of  a  series  hybrid  electric  vehicle  is  given  to  validate  the  genetic  algorithm 
results. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Hybrid  Electric  Vehicles  (HEVs)  have  great  capabilities  as  new 
alternative  means  of  transportation.  The  advantages  of  HEVs  over 
the  conventional  vehicles  are  mainly  improved  fuel  economy  and 
reduced  emissions.  In  a  series  HEV,  traction  force  is  provided  only 
by  an  electric  motor.  The  driving  power  for  this  motor  comes  from 
a  bidirectional  storage  system  (e.g.,  battery,  flywheel,  or  ultra¬ 
capacitor)  and  an  engine-generator  set  (Genset).  Now,  the  question 
is  how  to  distribute  the  required  power  for  the  electric  motor 
between  the  batteries  and  the  Genset.  The  answer  is  provided  by 
the  control  strategy  that  is  being  adopted.  A  control  strategy  is  an 
algorithm  regulating  the  operation  of  the  drive-train.  It  takes  data 
from  the  vehicle  (e.g.,  speed,  acceleration  and  grade)  and  makes 
decisions  to  turn  on/off  certain  components  or  to  increase/decrease 
their  power  output  [1  ]. 

Not  many  approaches  exist  in  literature  for  the  control  strat¬ 
egy  of  series  HEVs.  Different  methods  of  optimal  control  have  been 
used  in  these  papers  to  minimize  the  fuel  consumption  or  emis¬ 
sion.  Among  them,  we  can  refer  to  Refs.  [2,3].  Other  studies,  e.g., 
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Refs.  [4,5],  have  considered  the  internal  battery  losses  empirically 
or  as  an  average  value. 

Nevertheless,  the  batteries  performance  and  life  length  have 
not  been  addressed  in  most  proposed  algorithms.  In  order  to 
reduce  the  costs  associated  with  loss  of  life  cycle  in  batteries, 
their  charge/discharge  patterns  have  to  be  managed.  For  this  pur¬ 
pose,  it  has  been  recommended  to  avoid  deep  discharging  and 
frequent  charging,  since  this  will  cause  the  batteries  life  to  dete¬ 
riorate  dramatically.  The  main  objective  of  this  paper  is  to  propose 
an  algorithm  by  which  the  demanded  power  (from  the  main  vehicle 
controller)  gets  distributed  between  the  batteries  and  the  Genset 
in  such  a  way  that  the  batteries’  life  losses  are  as  low  as  possi¬ 
ble.  Besides,  in  the  proposed  control  strategy,  the  factors  reducing 
the  battery  life  such  as  deep  discharge  and  frequent  charging  are 
avoided.  This  looks  like  solving  an  optimization  problem  for  min¬ 
imizing  of  the  fuel  consumption  with  the  battery  life  expectation 
constraint.  The  method  will  be  verified  by  simulation. 

The  vehicle  under-study  which  is  used  for  current  investigations 
is  the  0457  bus  which  has  been  converted  to  a  series  hybrid  vehi¬ 
cle  at  Vehicle,  Fuel  and  Environment  (VFE)  Research  Institute  at 
University  of  Tehran. 

The  proposed  algorithm  relies  on  parameters  that  characterize 
the  driving  schedule  and  is  also  flexible  to  unexpected  changes 
in  the  driving  situation.  The  layout  of  this  paper  is  as  follows.  In 
Section  2,  the  series  HEV  model  and  power  losses  due  to  each  com- 
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Nomenclature 

SoC  battery  state  of  charge 

DoD  battery  depth  of  discharge 

R[  battery  internal  resistance  (£2) 

/b  battery  current  (A) 

Pb  battery  power  (W) 

Ub  battery  voltage  (V) 

U0  initial  capacitor  voltage  (V) 

Pd  drive  power  (W) 

Pe/g  Genset  power  (W) 

P*ls  desired  Genset  power  (W) 

a0,  a\ ,  a2  parameters  of  Genset  efficiency  curve 
Dr  DOD  for  which  rated  cycle  life  was  determined 
CR  rated  ampere-hours  capacity  at  rated  discharge  cur¬ 

rent  /R  (Ah) 

Lr  cycle  life  at  rated  DOD,  DR  and  discharge  current  /r 

(Ah) 

Da,  Ca  actual  value  of  DR  and  Cr 

u0,  U\ ,  u2  parameters  of  cell  cycle  life  curve 

Ltime  life  time  (s) 

^actual  actual  absolute  discharge 

de ff  effective  discharge 

Peff  effective  discharge  power  (W) 

Pb/ioss  power  losses  of  battery  ( W) 

Pbus  bus  power  (W) 

Hsoc  higher  value  of  SoC 

Lsoc  lower  value  of  SoC 

tmoff  shortest  allowed  duration  of  engine  off  time 

Poff  power  below  which  engine  is  shut  off  if  the  batteries 
are  charged 

Greek  symbols 

rjgen  efficiency  of  generator 

tjice  efficiency  of  ICE 

rR  rated  charge  life  (Ah) 


ponent  have  been  introduced.  The  optimization  problem  definition 
and  control  algorithm  have  been  explained  in  Section  3.  In  Section 
4,  results  from  the  developed  control  strategy  have  been  com¬ 
pared  to  the  common  control  strategy.  Conclusions  are  presented  in 
Section  5. 


2.  System  components 

The  series  HEV  power-train  is  composed  of  an  engine/generator 
set  (Genset),  an  electric  traction  motor  with  a  power  converter,  and 
a  battery  pack  with  an  appropriate  controller.  In  addition,  there  is  a 
central  controller  that  manages  the  interaction  between  these  com¬ 
ponents,  as  shown  in  Fig.  1 .  The  output  of  the  generator  is  connected 
to  an  electric  power  bus  through  a  Power  Conditioner  (PC).  The  bat¬ 
tery  pack  serves  as  the  bidirectional  electrochemical  energy  source 
which  is  connected  to  the  bus  by  means  of  a  power  electronic  con¬ 
verter  (DC/DC  converter).  The  electric  power  bus  is  also  connected 
to  the  controller  of  the  electric  traction  motor.  The  traction  motor 
can  be  controlled  either  as  a  motor  or  a  generator,  and  either  in 
forward  or  reverse  motion. 

To  study  the  optimization  algorithm  by  which  the  parameters 
of  the  control  strategy  are  determined,  the  model  of  a  series  hybrid 
electric  vehicle  was  built.  This  model  enables  us  to  calculate  the  flow 
of  power  among  different  components  of  the  drive-train.  Accuracy 
of  the  model  is  of  course  the  main  concern,  but  as  long  as  consistent 
behavior  is  seen,  a  few  percents  of  error  in  the  fuel  consumption  or 
battery  losses  are  tolerated. 

2.1.  Genset 

The  Genset  provides  one  part  of  the  power  that  is  required  by  the 
electric  traction  motor  for  propulsion.  As  proposed  in  Ref.  [4],  the 
convex  second  degree  polynomial  of  Genset’s  power  can  represent 
the  efficiency  curves  at  different  speeds,  as  follows  (Fig.  2): 

*7lCE  =  Oo  +  <3l  pe/g  +  a2P\/g  ( 1 ) 

where  Pe/g  is  the  mechanical  power  on  the  Genset  output  shaft 
and  r] ice  denotes  the  Internal  Combustion  Engine  (ICE)  thermal 
efficiency.  The  generator  efficiency  map  is  depicted  in  Fig.  3.  The 
reference  engine  speed  is  obtained  from  the  efficiency/power  map 
to  optimize  the  efficiency  at  the  requested  power.  Then  the  ICE 
controller  regulates  the  engine  shaft  speed  to  follow  the  reference 
value,  as  shown  in  Fig.  4. 

2.2.  Battery  pack 

The  battery  pack  is  an  energy  source  just  like  a  fuel  tank.  How¬ 
ever,  its  high  cost,  low  energy  storage  density  and  shorter  life  often 


Genset 


Gen:  Generator  PC:  Power  Conditioner  EM:  Electric  Motor 


Fig.  1.  General  arrangement  of  a  series  HEV. 
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Qmin(minimum  speed) 


Q=ICE  Speed 


Envelope  of  curves 


Fig.  2.  Qualitative  efficiency  of  the  ICE.  The  specific  consumption  curves,  obtained  from  the  engine  efficiency  maps  have  the  qualitative  behavior  of  the  left  graph,  one  curve 
for  each  rotational  speed.  To  ensure  the  best  efficiency  for  each  value  of  the  desired  power  of  Genset,  the  ICE  should  vary  its  speed  so  as  to  follow  the  envelope  of  such  curves, 
as  shown  in  the  right  graph. 


Fig.  3.  Efficiency  map  of  the  generator. 


makes  the  battery  pack  the  weakest  chain  ring  in  the  drive-train 

[6]. 

Batteries  are  expensive  items.  The  cost  for  a  certain  storage 
capacity  can  vary  from  100  to  900$kW_1  h-1  or  even  higher  [7]. 
Unlike  the  fuel  tank,  batteries  might  not  last  for  the  whole  life  span 
of  the  vehicle  and  their  condition  usually  deteriorates  over  their  life 
as  well.  The  life  of  a  battery  is  measured  in  cycles  or  years  for  cyclic 
applications  (like  in  HEVs).  The  cycle  life  of  batteries  is  between  300 
and  1500  full  cycles,  depending  on  the  battery  type. 

In  order  to  prolong  batteries’  life,  different  factors  of  their  oper¬ 
ation  can  be  taken  into  account.  However,  the  first  step  towards 
achieving  this  goal  is  to  implement  a  proper  life  prediction  method 
for  them.  The  method  used  in  this  paper  is  similar  to  the  method 
reported  in  Ref.  [6].  Each  battery  cell  has  a  finite  life  as  measured  by 


Fig.  4.  Control  arrangement  to  adjust  the  shaft  speed  to  the  requested  output  power. 
The  ICE  should  operate  so  that  it  follows  the  minimum  consumption  curve  by  varying 
its  rotational  speed  between  the  minimum  and  maximum  allowed  when  following 
the  required  power  P* 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Depth  of  Discharge  (DoD) 


Fig.  5.  Best  fit  curve  to  battery  manufacturer’s  cycle  life  data. 

the  sum  of  the  effective  ampere-hour  during  its  useful  life.  When 
the  cumulative  effective  ampere-hour  (the  total  of  the  individ¬ 
ual  effective  ampere-hour  corresponding  to  a  series  of  discharge 
events)  equal  the  rated  charge  life  of  the  cell,  the  cell  will  reach  its 
useful  life.  The  rated  charge  life,  r r,  of  a  cell  is  defined  as: 

rR  =  LrDrCr  (2) 

where,  Cr  is  the  rated  Ampere-hour  capacity  at  rated  discharge  cur¬ 
rent  /R.  Besides,  DR  stands  for  the  Depth  of  Discharge  (DoD)  for 
which  rated  cycle  life  was  determined,  and  LR  represents  the  cycle 
life  at  rated  DoD,  DR  and  discharge  current  /R. 

The  actual  charge  life  (cumulative  effective  ampere-hour)  of  a 
cell,  is  a  function  of  its  DoD.  To  determine  this  relationship,  the 


Current  (A) 

Fig.  6.  Actual  capacity  versus  discharge  current  for  a  68  Ah  valve  regulated  lead  acid 
battery. 
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R/SoCJ) 


Table  2 

Description  of  control  strategy  parameters  used  in  optimization. 


Optimization  variable 

Description 

Hsoc 

Highest  desired  battery  SoC 

Lsoc 

Lowest  desired  battery  SoC 

Po  ff 

Power  below  which  engine  is  shut  off  if  the  batteries 
are  charged 

Pch 

Power  applied  to  engine  to  charge  the  batteries 

hnoff 

Shortest  allowed  duration  of  engine  off  time 

Therefore,  the  total  power  losses  of  the  battery  are  obtained  as: 


Fig.  7.  Battery  equivalent  circuit  utilized  in  the  paper  for  modeling  of  the  energy 
storage  system. 


following  equation  is  proposed: 

L  =  U2^y°e»i(i-(D/DR))  (3) 

Fig.  5,  shows  the  best  fitted  curve  to  a  set  of  four  data  points: 
cycle  life  versus  DoD  which  can  be  provided  by  the  manufacturer 
of  a  lead-acid  battery,  based  on  Eq.  (3).  For  the  provided  battery,  the 
parameters  in  Eq.  (3)  are  obtained  as  u0  =  0.1,  U\  =  1.69  and  u2  =  765. 

It  is  assumed  that  for  a  given  actual  discharge,  the  effective  dis¬ 
charge  de ff,  will  increase  with  discharge  rate  and  can  be  roughly 
expressed  by  the  following  function: 


deff  = 


Cr 

Ca 


^actual 


(4) 


in  which,  CA  can  be  obtained  easily  from  the  manufacturer’s  data 
such  as  the  one  shown  in  Fig.  6. 

The  effects  of  discharge  rates  are  combined  simply  by  multiply¬ 
ing  the  factors  expressed  in  Eqs.  (3)  and  (4). 

de ff=  (^)U°eu'((DA/DR)-1)^dactual  (5) 


^b/loss  —  R'b+Peff  =  (l+M)R/2  (8) 

where, 

"  l)  °e-Ul(CA/cR)  (9) 

and,  RI^  is  the  losses  caused  by  internal  resistance  and  Peff  is  the 
battery  losses  reducing  the  battery  life. 

A  dynamical  model  for  the  lead-acid  battery  should  be  obtained 
for  the  purpose  of  minimization  [8].  Battery  modeling  is  of  major 
concern  in  the  overall  HEV  model  development.  Studies  have  been 
made  to  model  the  power  output  for  batteries  in  various  approaches 
[9,10].  A  simple  model  represented  in  Fig.  7,  as  an  equivalent  elec¬ 
trical  circuit,  is  considered  for  modeling  of  the  battery.  Battery 
parameters  which  are  used  for  modeling  are  depicted  in  Table  1. 
This  model  has  successfully  been  used  for  lead-acid  batteries  [9]. 

It  should  be  noted  that  both  the  internal  resistance  and  the 
electromotive  force  are  not  constants  during  the  charge/discharge 
processes  and  vary  as  a  function  of  the  State  of  Charge  (SoC)  and 
electrolyte  temperature.  For  simplicity,  these  dependencies  are 
assumed  to  be  linear. 

3.  Evolutionary  design  of  the  control  strategy 


It  should  be  noted  that  Eq.  (5)  expresses  the  effective  discharge  for  a 
single  discharge  of  specific  magnitude  and  rate.  A  life  prediction  for 
a  cell  subjected  to  an  irregular  pattern  of  charge-discharge  cycles 
requires  the  summation  of  the  effective  discharges  from  a  series  of 
discharge  events. 

The  discharge  profile  can  be  obtained  by  actually  monitoring  of 
the  battery  current  on  an  operating  system,  or  by  modeling  battery’s 
usage  on  a  proposed  system.  The  prescribed  series  of  n  discharge 
events  will  correspond  to  a  certain  time  period  (T)  of  system’s  oper¬ 
ation.  The  life  time  (Ltime)  °f  the  cell  under  the  specified  usage 
pattern  is  then  given  by: 


t'time  — 


A  _  LRDRCR 
reff/T  y^=  i^eff 


(6) 


In  order  to  formulate  the  battery  losses  that  reduce  its  life,  the 
following  equation  can  be  employed: 


Table  1 

Relevant  battery  data  used  in  simulation  of  the  HEV. 


Battery  rated  capacity  (Cr) 

68  Ah 

Battery  internal  resistance  R 

0.5  m£2 

Rated  DOD 

0.8 

Number  of  battery  module 

54 

Cycle  life  curve:  u o 

0.19 

Cycle  life  curve:  u\ 

1.69 

Cycle  life  curve:  U2 

765 

Efficiency  curve:  ao 

0.3520 

Efficiency  curve:  a\ 

0.0025  kW”1 

Efficiency  curve:  02 

-0.0001  kW“2 

Genetic  algorithm  (GA)  used  in  this  paper  as  the  evolutionary 
algorithm  is  motivated  by  the  mechanism  of  natural  selection,  a 
biological  process  in  which  fitter  individuals  have  higher  ability  to 
survive  during  evolution. 

The  principles  of  multi-objective  optimization  are  different 
from  that  of  a  single-objective  optimization.  The  main  goal  in  a 
single-objective  optimization  is  to  find  the  global  optimal  solu¬ 
tion.  However,  in  a  multi-objective  optimization  problem,  there 
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Table  3 

Relevant  vehicle  and  component  data  used  in  series  hybrid  electric  vehicle  model. 


Vehicle 

Vehicle  mass:  9100  kg 

27  x  70  kg  (according  to  SAE  931788) 
Front  area:  7.23  m2 

Rolling  resistance  coef.:  0.008 
Aerodynamic  coef.:  0.79 

Engine 

6.8  L  Cl,  170  l<W/2200  rpm 

Electric  motor 

Rated  power:  170  kW 

Max.  power:  300  kW 

Max.  torque:  1000  Nm  at  500  A 

Peak  efficiency:  94.1% 

PM  Generator 

Rated  power:  105  kW 

Max.  power:  160  kW 

Max.  torque:  450  Nm  at  250  A 

Peak  efficiency:  95.2% 

is  more  than  one  objective  function,  each  of  which  may  have  a 
different  individual  optimal  solution.  If  there  is  sufficient  differ¬ 
ence  in  the  optimal  solutions  corresponding  to  different  objectives, 
the  objective  functions  are  often  known  as  being  conflicting 
to  each  other.  Multi-objective  optimization  with  such  conflict¬ 
ing  objective  functions  gives  rise  to  a  set  of  optimal  solutions, 
instead  of  one  optimal  solution.  The  reason  for  the  optimal¬ 
ity  of  many  solutions  is  that  no  solution  can  be  considered  to 
be  better  than  the  other  with  respect  to  all  objective  functions. 
These  optimal  solutions  have  a  special  name— Pareto-optimal  solu¬ 
tions. 

It  is  clear  that  the  concept  of  optimality  in  multi-objective  opti¬ 
mization  deals  with  a  number  (or  a  set)  of  solutions,  instead  of  one 
solution.  Based  on  the  above  discussions,  the  conditions  for  a  solu¬ 
tion  to  become  dominated  with  respect  to  another  solution  and 
then  present  conditions  for  a  set  of  solutions  to  become  a  Pareto- 
optimal  set  are  defined. 

Let  us  assume  that  the  set  of  variables  x*,  i  =  1,. . .,  M,  is  a  set 
of  decision  variables  for  a  multi-objective  optimization  problem.  In 
order  to  measure  the  decency  of  a  certain  solution,  some  criteria  for 
evaluating  the  solution  quality  should  be  defined.  These  criteria  are 
expressed  as  a  set  of  functions/i  (x),. .  .,/j(x),. .  .,/p(x)  of  the  decision 
variables,  which  are  called  objective  functions.  Oftentimes,  some 
of  them  are  contradicting  with  others.  The  constraints  specify  the 
feasible  region  X  in  which  any  point  x  defines  a  feasible  solution. 
It  is  impossible  that  all  of  the  /j(x),  j  =  1,. . .,  P  values  have  an  opti¬ 
mum  in X  at  a  common  point  x.  Certain  criteria  should  be  developed 
to  determine  an  optimal  solution  in  this  circumstance.  One  inter¬ 
pretation  of  the  term  optimum  in  the  multi-objective  optimization 
scenarios  is  the  Pareto-optimality,  which  are  highly  associated  with 
the  concept  of  dominance.  A  solution  Xi  dominates  a  solution  x2  if 
and  only  if  the  two  following  conditions  are  satisfied  at  the  same 


Fig.  9.  Illustration  of  thermostat  control  strategy.  When  the  SoC  of  the  battery 
reaches  its  preset  top  line,  the  Genset  is  turned  off  and  will  be  turned  on  when 
the  SoC  of  the  battery  reaches  its  bottom  line. 


Table  4 

Comparison  of  main  parameters  obtained  from  optimization  and  equivalent  value 
of  them  in  thermostatic  control  strategy. 


Parameters 

Developed  control  strategy 

Thermostatic 

HSoc 

0.915 

0.8 

f-SoC 

0.35 

0.3 

A>ff 

35  kW 

75  kW 

Pch 

28  kW 

0 

bnoff 

355 

00 

Table  5 

Minimum  requirements  of  road  performance  based  on  previous  version  of 0457  city. 

Max  vehicle  speed  (km  h-1 )  110 

Max  vehicle  speed  at  6%  grade  (km  h- 1 )  85 

Time  to  accelerate  from  0  to  100  km  h_1  (s)  35 


time:  firstly,  X\  is  no  worse  than  x2  in  all  objective  evaluations,  i.e., 
fj{x\ )  <fj(x2)  for  all  j  =  1,. . .,  P.  Secondly,  Xi  is  strictly  better  than  x2 
in  at  least  one  objective,  i.e.,/j(xi  )</j(x2)  for  at  least  one  j  e  {1,. . ., 
P}.  As  a  result,  in  a  set  of  Pareto-optimal  solutions,  there  is  no 
solution  which  dominates  another  with  respect  to  all  the  design 
objectives  involved.  It  should  be  noted  that  multi-objective  opti¬ 
mization  needs  a  decision-making  process  as  there  is  not  a  single 
solution  but  a  set  of  non-dominated  solutions  from  which  the  best 
should  be  chosen.  That  is,  the  major  two  tasks  of  multi-objective 
optimization  are  to  obtain  a  representative  set  of  non-dominated 
solutions  and  then  select  a  suitable  solution  from  this  set  based  on 
the  specific  criterion.  In  this  study,  the  standard  GA  algorithm  is 
improved  and  extended  to  handle  the  target  application,  which  is 
essentially  a  multi-variable  problem. 

3.1.  Problem  formulation 

Optimization  objective  is  to  reduce  the  losses  of  the  system, 
which  are  reflected  by  four  parameters  including  the  power  losses 
of  ICE,  generator,  battery  pack  and  the  losses  due  to  the  battery  life. 
According  to  the  above  formulation,  the  set  of  objective  functions 
are  chosen  as: 


/i(x)  =  (§)'  U°  (§  -  i)  VUl(CA/CR)R/g 

(10) 

(ID 

A  00  =  (1  —  *?gen  )^e/g 

(12) 

f4(x)  =  R/2 

(13) 

It  should  be  mentioned  that  due  to  the  direct  dependence  of  the 
electric  traction  motor’s  torque  and  speed  on  the  driver  commands 
and  road  conditions,  the  power  losses  due  to  the  components  can 
be  ignored. 

3.2.  Archiving 

The  major  function  of  the  archive  is  to  store  a  historical  record  of 
the  non-dominated  solutions  found  along  the  heuristic  search  pro¬ 
cess.  The  archive  interacts  with  the  generational  population  at  each 
iteration  so  as  to  absorb  superior  current  non-dominated  solutions 


Table  6 

Comparison  of  the  results  for  the  developed  and  thermostatic  control  strategy. 


Fuel  consumption 

Battery  energy 

Battery  life  time 

(1100  km-1) 

losses  (kj) 

(no.  of  cycles)3 

Developed  controller 

65.4 

897 

1843 

Thermostat  control 

71.5 

1352 

1360 

a  This  means  that  how  many  cycles  batteries  can  work  before  failure  of  health. 


M.  Amiri  et  al.  /  Journal  of  Power  Sources  190  (2009)  372-379 


377 


Fig.  10.  Tehran  city  bus  driving  cycle  test  results  using  developed  control  strategy,  (a)  Efficiency  map  of  the  electric  motor,  (b)  Electric  motor  output  power,  (c)  ICE  output 
power,  (d)  ICE  efficiency  history,  (e)  Battery  power  loss,  (f)  SoC  history. 


and  eliminate  inferior  solutions  currently  stored  in  the  archive.  The 
non-dominated  solutions  obtained  at  every  iteration  in  the  gener¬ 
ational  population  are  compared  with  the  contents  of  archive  in  a 
one-per-one  basis.  A  candidate  solution  can  be  added  to  the  archive 
if  it  meets  any  of  the  following  four  conditions: 


•  There  is  no  solution  currently  stored  in  the  archive; 

•  The  archive  is  not  full  and  the  candidate  solution  is  not  dominated 
by  or  equal  to  any  solution  currently  stored; 

•  The  candidate  solution  dominates  any  existing  solution  in  the 
archive; 
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Fig.  11.  Tehran  city  bus  driving  cycle  test  results  based  on  thermostatic  controller,  (a)  ICE  output  power,  (b)  ICE  efficiency  history,  (c)  Battery  power  loss,  (d)  SoC  history. 


•  The  archive  is  full  but  the  candidate  solution  is  non-dominated 
and  is  in  a  sparser  region  than  at  least  one  solution  currently 
stored. 

During  the  heuristic  multi-objective  optimization  process,  the 
multi-objective  optimization  algorithm  is  attempting  to  build  up  a 
discrete  picture  of  a  possibly  continuous  Pareto  front.  Therefore  it  is 
often  desired  to  distribute  the  solutions  as  diversely  as  possible  on 
the  discovered  tradeoff  curve.  Furthermore,  the  uniformity  among 
the  distributed  solutions  is  also  crucial  so  as  to  achieve  consistent 
and  smooth  transition  among  the  solution  points  (when  searching 
for  the  preferred  solution  based  on  the  particular  requirements  of 
the  target  problem).  Therefore,  to  accomplish  these  challenges,  it 
is  necessary  to  preserve  the  diversity  of  the  solutions  distribution 
during  the  optimization  process. 

3.3.  Encoding  scheme 

It  is  significant  to  select  the  appropriate  optimal  variables  for 
the  optimization  process.  The  simplest  way  to  implement  a  control 
strategy  for  the  power  management  of  an  HEV  is  to  propose  a  set 
of  static  thresholds,  described  in  details  in  Table  2.  The  selected 
variables  correlate  closely  with  the  fundamental  vehicle  operation. 

In  adopting  GAs  to  the  scheduling  problem  such  as  one  dealt 
in  this  paper,  each  chromosome  represents  a  candidate  solution 


consisting  of  the  following  genes: 

(f^SoC?  f'SoC?  ^off>  ^ch>  hnoff)  (13) 

Each  element  in  the  chromosome  is  coded  using  a  floating-point 
number. 

4.  Performance  simulation 

In  order  to  help  with  the  simulation,  design,  and  analysis  of 
the  designed  power-train,  certain  typical  driving  schedules  have 
been  developed.  These  driving  schedules  represent  typical  traffic 
conditions  for  a  particular  range  of  time.  The  optimization  process 
is  operated  using  Tehran  transit  bus  driving  cycle  shown  in  Fig.  8 
[11  ].  The  total  test  time  for  this  cycle  is  1800  s,  the  average  speed  is 
9.61  kmh-1  and  the  maximum  speed  is  50.11  km  h-1.  The  distance 
driven  is  approximately  4.81  km. 

After  watchful  tuning,  the  settings  of  GA  parameters  in  the  sim¬ 
ulation  are  as  follows:  number  of  generations  is  432,  population 
size  is  170,  crossover  probability  is  87%,  and  mutation  probability 
is  1.5%. 

In  order  to  investigate  the  advantages,  the  developed  control 
strategy  was  implemented  in  Hybrid  Electric  Vehicle  FEed-foRwarD 
SIMulator  (HEV-FERDSIM),  developed  in  VFE  research  institute  [  12 ]. 
The  simulation  parameters  for  vehicle  used  and  its  components  are 
listed  in  Table  3.  The  battery  SoC  correction  procedure  is  used  to 
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correct  the  fuel  economy  for  the  case  in  which  the  initial  and  the 
final  battery  SoC  are  not  the  same  [12]. 

An  optimized  control  strategy  called  thermostatic  (engine- 
on-off)  is  used  for  comparison,  illustrated  in  Fig.  9.  The  operation  of 
the  engine/generator  is  totally  controlled  by  the  SoC  of  the  battery 
pack.  When  the  SoC  of  the  batteries  reaches  its  predetermined  high 
value,  the  engine/generator  is  turned  off  and  the  vehicle  is  driven 
only  by  the  battery  pack.  On  the  other  hand,  when  the  SoC  of  the  bat¬ 
teries  reaches  its  predetermined  low  value,  the  engine/generator  is 
turned  on.  The  battery  pack  charges  by  the  engine/generator.  Thus, 
the  engine  can  always  operate  within  its  optimal  region. 

The  control  strategy  parameters  obtained  through  the  GA  opti¬ 
mization  are  shown  in  Table  4,  in  addition  to  thermostatic  control 
strategy  parameters.  It  should  be  noted  that  some  design  param¬ 
eters  are  not  directly  derived  from  the  optimization,  but  are  cal¬ 
culated  from  other  design  parameters.  Each  control  strategy  must 
be  able  to  drive  the  mentioned  transit  bus,  in  such  a  way  that  the 
minimum  requirements  of  performance  (i.e.,  Table  5),  are  satisfied. 

Table  6  presents  fuel  consumption,  battery  energy  losses  and 
battery  life  time  for  the  developed  and  thermostatic  controllers.  As 
it  can  be  seen,  not  only  the  fuel  consumption  is  decreased  when 
using  the  developed  controller,  but  also  the  battery  energy  losses 
and  life  are  improved  considerably. 

Simulation  results  of  the  vehicle  under  the  prescribed  controller 
policy  are  shown  in  Fig.  10.  The  electric  motor  operating  points 
and  power  output  are  presented  in  Fig.  10(a)  and  (b),  respectively. 
The  ICE  power  output  and  its  efficiency  points  during  the  cycle 
are  depicted  in  Fig.  10(c)  and  (d),  respectively.  As  can  be  seen, 
the  operating  points  in  Fig.  10(d)  are  close  to  the  optimal  region, 
which  indicates  that  the  ICE  has  been  operating  close  to  maximum 
possible  efficiency. 

The  battery  power  losses  and  SoC  history  during  the  cycle  are 
shown  in  Fig.  10(e)  and  (f).  It  can  be  concluded  that  lower  variation 
of  battery  SoC  during  the  cycle  is  the  main  reason  of  increasing  of 
the  battery  life. 

In  the  case  of  using  thermostatic  controller  in  the  simulation,  the 
battery  power  losses,  SoC  history,  power  output  and  efficiency  are 
shown  in  Fig.  11  (a-d),  respectively.  As  can  be  seen,  the  ICE  is  turned 
off  most  of  the  time  and  battery  is  fully  discharged.  Therefore,  the 
battery  power  losses  increase  in  comparison  to  the  developed  con¬ 
troller  losses. 

5.  Conclusions 

A  general  control  strategy  for  series  FIEVs  has  been  proposed  in 
order  to  optimize  power-train  efficiency  based  on  the  minimiza¬ 


tion  of  the  losses  reducing  battery  life.  As  long  as  ICE  efficiency 
is  high  enough,  there  will  be  no  need  for  charging  or  discharging 
the  batteries  which  is  an  inefficient  process.  This  control  strategy 
has  been  evaluated  in  detail  by  means  of  validated  simulations 
of  an  FIEV.  Since  the  simulations  have  been  performed  based 
on  an  existing  vehicle,  the  intention  of  the  authors  is  to  make 
experimental  validations  of  the  proposed  techniques  in  the  next 
step.  Because  of  its  great  importance  in  the  vehicle  control  and 
its  inherent  complexity,  more  study  is  required  to  address  chal¬ 
lenging  technical  issues  of  the  algorithm  and  the  whole  control 
strategy. 
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